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ABSTRACT 

We present 230 GHz Submillimeter Array continuum and molecular line observations of the newly 
discovered FUor candidate HBC722. We report the detection of seven 1.3 mm continuum sources 
in the vicinity of HBC722, none of which correspond to HBC722 itself. We compile infrared and 
submillimeter continuum photometry of each source from previous studies and conclude that three are 
Class embedded protostars, one is a Class I embedded protostar, one is a Class I/II transition object, 
and two are either starless cores or very young, very low luminosity protostars or first hydrostatic 
cores. We detect a northwest-southeast outflow, consistent with the previous detection of such an 
outflow in low-resolution, single-dish observations, and note that its axis may be precessing. We 
show that this outflow is centered on and driven by one of the nearby Class sources rather than 
HBC722, and find no conclusive evidence that HBC722 itself is driving an outflow. The non-detection 
of HBC722 in the 1.3 mm continuum observations suggests an upper limit of 0.02 M for the mass 
of the circumstellar disk. This limit is consistent with typical T Tauri disks and with a disk that 
provides sufficient mass to power the burst. 

Subject headings: stars: formation - stars: low-mass - stars: protostars - stars: flare - stars: variables: 
T Tauri, Herbig Ae/Be - ISM: individual objects (HBC722) - ISM: jets and outflows 



1. INTRODUCTION 

FU Orionis objects (hereafter FUors) are a group 
of young, pre-main sequence stars observed to flare in 
brightness by 4 — 6 magnitudes in the optical and re- 
main bright for decades (Herbig 1977). They are named 
after the prototype FU Orionis, which flared by about 
6 magnitudes in 1936 and has remained in an elevated 
state to the present day (Wachmann 1954; Herbig 1966). 
Only 10 confirmed FUors are known to exist from di- 
rect observations of flares, with about another 10 identi- 
fied based on similar spectral characteristics to the con- 
firmed FUors (see Reipurth & Aspin 2010 for a recent 
review). The large amplitude flares are attributed to 
enhanced accretion from the surrounding circumstellar 
disk (Hartmann & Kenyon 1985), with the accretion rate 
from the disk onto the star increasing to up to ~ 10 -4 
M yr _1 (Hartmann & Kenyon 1996). Various trig- 
gering mechanisms for the accretion bursts have been 
proposed, including interactions with binary companions 
(Bonnell & Bastien 1992), thermal instabilities (Hart- 
mann & Kenyon 1996), and gravitational and magne- 
torotational instabilities (Zhu et al. 2007, 2009a, 2009b; 
Vorobyov & Basu 2010). FUors are especially interest- 
ing and relevant to general star formation studies because 
they may represent the late, optically visible end stages 
of episodic accretion bursts and luminosity flares through 
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the duration of the embedded phase (e.g., Kenyon et 
al. 1990; Enoch et al. 2009; Evans et al. 2009; Vorobyov 
2009; Dunham et al. 2010; Dunham & Vorobyov 2012). 
With so few bona-fide FUors, detailed observations of 
each one are necessary in order to characterize their prop- 
erties and understand their place in the general star for- 
mation process. 

In this paper, we present 230 GHz Submillimeter Array 
(SMA; Ho et al. 2004) observations of the newly discov- 
ered FUor candidate HBC722. As described in more de- 
tail in $2] below, HBC722 is located within a small group 
of ~ 10 young stars, greatly complicating the analysis 
of existing, low-resolution single-dish ground and space- 
based (sub) millimeter data. The SMA 230 GHz contin- 
uum and molecular line observations presented in this 
paper are motivated by a need to disentagle the millime- 
ter emission from the various sources in the vicinity of 
HBC722 in order to better determine its evolutionary 
status and physical properties. The organization of this 
paper is as follows: a brief summary of HBC722 is given 
in $2j a description of the observations and data reduc- 
tion is provided in £[3j the basic results are presented in 
includin g th e continuum data in §4.11 and the CO 
line data in §4.2[ a dis cussi on of the detected continuum 
sources is presented in §5.1[ a discussion of the evolution- 
ary status of HBC722 is given in l5.2[ and a summary of 
our results is presented in £0 

2. HBC722 

HBC722, also known as V2493 Cyg, LkHa 188-G4, 
and PTFlOqpf, is located at RA. = 20:58:17.03, decl. = 
+43:53:43.4 (J2000) in the "Gulf of Mexico" region of the 
North American/Pelican Nebula Complex at a distance 
of 520 pc (Straizys et al. 1989; Laugalys et al. 2006). 
Prior to 2010 it was regarded as an emission- line star 
with a spectral type of K7— M0 and Ay of 3.4 magni- 
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tildes (Cohen & Kuhi 1979). Semkov et al. (2010) and 
Miller et al. (2011) independently reported a 4 — 5 mag- 
nitude optical flare in HBC722. This flare began some- 
time before 2010 May, reached peak brightness in late 
2010 September, followed a similar rise in brightness as 
other FUors, and exhibited spectral characteristics in- 
dicative of FUors (Semkov & Peneva 2010a; Semkov & 
Peneva 2010b; Semkov et al. 2010; Munari et al. 2010; 
Leoni et al. 2010; Miller et al. 2011). Prior to the flare, 
HBC722 featured a spectral energy distribution (SED) 
consistent with a Class II T Tauri star, with L^i = 0.85 
L and infrared spectral slope a = —0.77 (Kospal et 
al. 2011; Miller et al. 2011). During the flare L boi in- 
creased to 8.7 — 12 L . After the peak brightness was 
reached in late 2010 September, HBC722 entered a phase 
of rapid decline, decreasing by ~ 0.5 magnitudes in the 
optical by 2010 December (Kospal et al. 2011). Based 
on a linear extrapolation of this initial decline, Kospal et 
al. (2011) predicted a return to quiescence by late 2011 
or early 2012, much too quickly to be a FUor. They also 
noted that the outburst L^i and mass accretion rate 
implied by this L^i are on the low end for FUors. How- 
ever, unpublished photometry from the American Asso- 
ciation of Variable Star Observers (AAVSO) 7 shows that 
the optical brightness of HBC722 stopped its decline in 
early 2011, remained constant throughout most of 2011 
at 3 — 4 magnitudes above the quiescent brightness, and 
increased by about 0.5 — 1 magnitudes between 2011 
September and 2012 June. The cessation in early 2011 of 
the initial rapid decline has been very recently confirmed 
by Lorenzetti et al. (2012) and Semkov et al. (2012). 

HBC722 was imaged by the Spitzer Space Telescope 
(Werner et al. 2004) at 3.6 — 8 jam with the Infrared 
Array Camera (IRAC; Fazio et al. 2004) and 24 - 160 
jam with the Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004) as part of a large survey of the 
North American and Pelican Nebula Complex (Guieu et 
al. 2009; Rebull et al. 2011). HBC722 is located within a 
small group of ~ 10 young stars called the LkHa group by 
Cohen & Kuhi (1979), all located within approximately 
20 - 30" (10400 - 15600 AU at the assumed distance of 
520 pc) of HBC722. Several are detected in the mid- 
infrared with Spitzer and are thus still associated with 
dust in surrounding disks and envelopes, emphasizing the 
complex nature of this region (see Figure 1 of Green et 
al. [2011] and g5ZQ below). 

Green et al. (2011) presented (sub) millimeter data 
on HBC722 and its surrounding region observed with 
the Herschel Space Observatory (Pilbratt et al. 2011) 
and Caltech Submillimeter Observatory (CSO), includ- 
ing Herschel 70 — 500 jam images with 5 — 36" resolution, 
Herschel 50 — 600 jam full spectral scans, a 350 /am image 
with 9" resolution obtained with the Submillimeter High 
Angular Resolution Camera II (SHARC-II) at the CSO, 
and a 12 CO J = 2 — 1 map with 30" resolution obtained 
at the CSO. HBC722 is detected at 70 /im with Herschel; 
at all longer wavelengths the resolution is insufficient 
to resolve HBC722 from other, nearby sources. Green 
et al. noted that the 100 — 500 /im Herschel emission 
does not peak on HBC722, suggesting it is not the dom- 
inant source of submillimeter emission. No emission is 
detected at the position of HBC722 in higher-resolution 

7 Available at http://www.aavso.org/ 



350 jam SHARC-II data, suggesting that this source is 
no longer associated with a circumstellar envelope, con- 
sistent with the pre-outburst classification as a Class II 
T Tauri star. The 30" resolution 12 CO J = 2 - 1 map 
presented by Green et al. clearly shows a NW-SE out- 
flow centered near HBC722 (see their Figure 2 and §4.21 
below), but the resolution is insufficient to definitively 
identify the driving source. In this paper, we present 
high angular resolution 1.3 mm continuum and 12 CO J 
= 2 — 1 observations obtained with the SMA in order to 
identify which sources in the vicinity of HBC722 are as- 
sociated with (sub) millimeter continuum emission, char- 
acterize each source, and identify the driving source (s) 
of the outflow(s) in this region. 

3. OBSERVATIONS AND DATA REDUCTION 

One track of observations of HBC722 was obtained 
with the SMA on 20 May 2011 in the compact config- 
uration with seven antennas, providing projected base- 
lines ranging from 5 to 76 m. A two pointing mo- 
saic was adopted to fully map the multiple sources in 
the vicinity of HBC722 (see, e.g., Figure 1 of Green 
et al. [2011]). The phase centers of the two point- 
ings are R.A.=20:58:16.56, decl.=+43:53:52.9 (J2000) 
and R.A.=20:58:17.67, decl.=+43:53:31.0 (J2000), with 
the total observing time divided equally between the 
two pointings. The observations were obtained with the 
230 GHz receiver and included 4 GHz of bandwidth per 
sideband, with 12 GHz spacing between the centers of 
the two sidebands. The correlator was configured such 
that the lower sideband (LSB) covered approximately 
216.8 - 220.8 GHz while the upper sideband (USB) cov- 
ered aproximately 228.8 — 232.8 GHz, providing simulta- 
neous observations of the 12 CO, 13 CO, C 18 J = 2 - 1, 
N 2 D+ J = 3 - 2, and SiO J = 5 - 4 lines. The lines 
were observed with either 256 ( 13 CO, SiO) or 512 ( 12 CO, 
C 18 0, N 2 D+) channels in the 104 MHz bands, providing 
channel separations of 0.53 and 0.26 km s _1 , respectively. 
The remaining bands were used to measure the 1.3 mm 
continuum with a total bandwidth of 5.49 GHz. 

The observations were obtained in moderate weather 
conditions, with the zenith opacity at 225 GHz ranging 
between 0.2 — 0.25 and the system temperature typically 
~ 250 K, ranging between 150 — 400 K depending on el- 
evation. Regular observations of the sources MWC349a 
and BLLAC were interspersed with those of HBC722 for 
gain calibration. Saturn and 3C279 were used for pass- 
band calibration, and Uranus was used for absolute flux 
calibration. We estimate a 20% uncertainty in the ab- 
solute flux calibration by comparing the measured fluxes 
of the calibrators from our calibrated data with those in 
the SMA calibrator database 8 for the same observation 
date. 

The data were inspected, flagged, and calibrated us- 
ing the MIR software package 9 and imaged, cleaned, and 
restored using the Multichannel Image Reconstruction, 
Image Analysis, and Display (MIRIAD) software package 
configured for the SMA 10 . Since our data are a mosaic of 
two pointings, we originally tried to use the MIRIAD 
clean task mossdi appropriate for mosaics. However, 

8 Available at http://smal.sma.hawaii.edu/callist/callist.html 

9 Available at https://www.cfa.harvard.edu/~cqi/mircook.html 

10 Available at http://www.cfa.harvard.edu/sma/miriad/ 
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TABLE 1 
SMA Observations of HBC722 



v Beam FWHM Beam PA a Bandwidth SV h la rms c 

Line (GHz) (Arcseconds) (Degrees) (GHz) (km s _1 ) (mJy beam -1 ) 

12 CO J = 2 - 1 230.53797 2.93 x 2.74 -50.5 0.082 0.5 126 

13 COJ = 2-l 220.39868 3.15 x 2.81 -46.0 0.082 05 103 

Continuum 225.44882 2.91 x 2.74 -52.1 5.494 ••• 1.65 

a Position angle of the long axis of the beam, measured east (counterclockwise) from north. 
b Width of each channel in km s -1 . 

c For lines, the mean of the la rms of the spectrum at each spatial position with the spectral resolution given 
in the previous column. For the continuum, the la rms of the continuum intensity. 



T 




I , , , I , , , I 

20 h 58 m 20 s 20 h 58 m 18 s 20 h 58 m 16 s 

RA (J2000.0) 



Fig. 1. UKIDSS K s band image of the vicinity of HBC722, 
taken from Green et al. (2011). The cross marks the position of 
HBC722 itself. Overplotted with solid black lines are SMA 1.3 
mm emission contours starting at 2a and increasing by 2a, where 
the la rms in the continuum image is 1.65 mJy beam -1 . Also 
plotted with dotted black lines are SMA 1.3 mm emission contours 
starting at —2a and decreasing by 2a. The synthesized beam size 
and shape of the SMA 1.3 mm continuum observations is shown 
by the black filled oval in the lower right. 

mossdi only allows for a single cleaning iteration, and 
since the region surrounding HBC722 shows complicated 
morphology, especially in the 12 CO data (see §H below), 
we found the result to be unacceptable since emission 
from the sidelobes of the dirty beam were very obviously 
present in our final image independent of the exact clean- 
ing parameters adopted. Thus, we instead cleaned and 
imaged each of the two pointings separately using the 
MIRIAD task clean, using an iterative cleaning process 
where we first cleaned only those regions showing clear 
emission in the dirty maps and then used those results 
as input models to further passes of clean applied over 
the full images. In the first iteration the regions var- 
ied from one channel to the next in the molecular line 
data since the emission morphology is different at dif- 
ferent velocities. Finally, we mosaiced together the final 
images from each of the two pointings. Both pointings 
were corrected for primary beam attenuation before mo- 
saicing. While all quantitative results in this paper are 
derived from the mosaics created after such correction, 
for display purposes the images and figures are created 
from the mosaics created without such corrections, un- 
less otherwise indicated. Imaging was performed with a 



robust uv weighting parameter of +1. The final maps 
were re-gridded onto 0.5" pixels. 

Only the continuum and 12 CO and 13 CO are detected 
and discussed in this paper. Table Q] lists, for each of 
these observations, the frequency of observation, the syn- 
thesized beam size and orientation, the total bandwidth, 
the channel separation (for the lines), and the measured 
la rms. For the continuum, the la rms is determined by 
calculating the standard deviation of all off-source pixels. 
For the spectra, the la rms is determined by calculating, 
for each pixel, the standard deviation of the intensity in 
each spectral channel outside of the velocity range of —20 
to 30 km s _1 , and then calculating the mean standard 
deviation over all pixels. 

4. RESULTS 
4.1. Continuum 

Figure [1] displays the UKIDSS K s band image of the 
vicinity of HBC722, taken from Green et al. (2011). 
Overlaid are contours showing the SMA 1.3 mm con- 
tinuum intensity. Seven individual millimeter sources 
are detected and labeled as MMS1 - MMS7 in order 
of decreasing declination. As evident from this Fig- 
ure, HBC722 itself is not detected at 1.3 mm. Asso- 
ciations between the seven millimeter sources detected 
by the SMA and known sources at infrared and submil- 
limeter wavelengths will be discussed in ^5.11 and the 
non-detection of HBC722 itself will be discussed in £15.21 

For all seven sources, we used the MIRIAD task imfit 
to fit an elliptical Guassian to each source. The proper- 
ties of each source as derived from these fits, including 
the peak position, peak flux density, total flux density, 
deconvolved source size, and deconvolved position angle, 
are reported in Table [2j The last two columns of Table 
[2] list the mass and density of each source derived from 
the SMA continuum detections. The mass is calculated 
as 

t5 v \±r>)K v 

where S v is the total flux density, B u (Td) is the Planck 
function at the isothermal dust temperature Tp, k v is 
the dust opacity, d = 520 pc, and the factor of 100 is the 
assumed gas-to-dust ratio. We adopt the dust opacities 
of Ossenkopf & Henning (1994) appropriate for thin ice 
mantles after 10 5 yr of coagulation at a gas density of 
10 6 cm -3 (OH5 dust), giving k v = 0.867 cm 2 gm _1 at 
the frequency of the continuum observations. To calcu- 
late the mass we assume To = 30 K. No uncertainties are 
given for the masses because these uncertainties are dom- 
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TABLE 2 

Elliptical Gaussian Fits to the 1.3 mm Continuum Detections 





Peak R.A. 


Peak Deck 


Peak Flux Density a 


Total Flux Density b 


Source Size c 


Source P.A. C 


M d 




n d 




Source 


(J2000) 


(J2000) 


(mJy beam -1 ) 


(mJy) 


n 


(degrees) 


(M ) 


(cm 


- 3 ) 


MMS1 


20:58:16.52 


+43:53:53.7 


27.4 ± 6.0 


38.3 ± 7.7 


2.0 x 1.6 


+26.6 


0.15 


5.3 


X 


10 7 


MMS2 


20:58:16.78 


+43:53:36.0 


7.7 ± 2.1 


55.1 ± 11.0 


8.3 x 6.1 


-70.4 


0.21 


1.2 


X 


10 6 


MMS3 


20:58:17.70 


+43:53:29.2 


18.6 ± 4.1 


66.1 ± 13.2 


4.9 x 4.3 


+31.8 


0.25 


5.3 


X 


10 6 


MMS4 


20:58:16.32 


+43:53:25.5 


33.2 ± 7.1 


60.5 ± 12.1 


3.3 x 1.8 


+30.2 


0.23 


3.2 


X 


10 7 


MMS5 


20:58:16.01 


+43:53:22.8 


25.1 ± 5.5 


62.7 ± 12.5 


5.9 x 1.0 


+55.6 


0.24 


3.4 


X 


10 7 


MMS6 


20:58:15.62 


+43:53:18.1 


15.5 ± 3.8 


44.3 ± 8.9 


6.5 x 1.6 


+89.4 


0.17 


1.0 


X 


10 7 


MMS7 


20:58:17.06 


+43:53:15.0 


11.0 ± 2.5 


20.3 ± 4.1 


3.1 x 2.1 


+1.2 


0.08 


1.0 


X 


10 7 



a Uncertainties include the statistical uncertainty returned by imfit and the 20% calibration uncertainy added in quadrature. 

k Uncertainties include only the 20% calibration uncertainty since imfit does not return a statistical uncertainty for this parameter. 

Deconvolved with the beam (see Tabled] for the beam size and shape). The source sizes for the weakest two sources (MMS2 and MMS7) are 
likely quite uncertain and are best treated as estimates only. 

d No uncertainties are given for the mass and mean number density since they are dominated by the uncertain assumptions for the dust 
temperature and opacity. See the text in £j4.1| for more details. 
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Fig. 2. SMA 12 CO J = 2-1 channel maps of the region surrounding HBC722 with a velocity resolution of 1 km s 1 , with the velocities 
marked in the top left of each panel. The central (0,0) position is that of HBC722 itself and is marked with a cross in each panel. The 
positions of the 7 sources detected in the SMA 1.3 mm continuum data are plotted as triangles and labeled in the first panel. The gray 
scale shows intensity in units of Jy beam -1 with a linear scaling ranging from -0.05 Jy beam -1 (white) to 4.0 Jy beam -1 (black). The 
solid contours start at 5cr and increase by 10cr, where the la rms is 0.094 Jy beam -1 in the 1 km s -1 channels. The lowest two contours 
are plotted in black; all others are plotted in white. The dotted contours start at — 5cr and decrease by 10cr. The synthesized beam size 
and shape of the 12 CO J = 2 — 1 observations is shown by the black filled oval in the lower right of each panel. 



inated by uncertain assumptions for the dust tempera- 
ture and opacity. For example, changing the assumed 
Td to 10 K would increase the masses by a factor of 4.5. 
Furthermore, the dust opacity at 1.3 mm can vary by fac- 
tors of ~ 2 — 4 depending on which dust opacity model is 
adopted (e.g., Shirley et al. 2005; 2011), directly leading 
to factors of ~ 2 — 4 variation in the mass. 

The mean number density of each source, n, is calcu- 
lated assuming spherical symmetry as 



M 



47r/im H r 6 eff 



(2) 



where M is the mass, r e ff is the effective radius 11 , run is 
the hydrogen mass, and \i is the mean molecular weight 
per free particle. We adopt \i = 2.37 for gas that is 71% 

11 The effective radius is defined as the geometric mean of the 
semimajor and semiminor axes. 
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Fig. 3. — High- velocity redshifted and blueshifted integrated 
12 CO J = 2 — 1 emission overlaid on the SMA 1.3 mm continuum 
image of the region surrounding HBC722. The central (0,0) posi- 
tion and the black cross mark the position of HBC722 itself. The 
grayscale is inverted and displayed in a linear stretch ranging from - 
1.65 mjy beam -1 (— la, white) to 28.05 mJy beam -1 (17a, black). 
The solid green contours show the continuum intensity starting at 
3a and increasing by 2a, where the la rms in the continuum image 
is 1.65 mJy beam -1 . The blue contours show blueshifted 12 CO J 
= 2 — 1 emission integrated from —5 to 1 km s -1 , while the red 
contours show redshifted emission integrated from 8 to 14 km s -1 . 
The solid blue and red contours start at 2.0 Jy beam -1 km s -1 
and increase by 2.0 Jy beam -1 km s -1 whereas the dotted blue 
and red contours start at —2.0 Jy beam -1 km s -1 and decrease 
by 2.0 Jy beam -1 km s -1 . The synthesized beam size and shape 
of the 1 CO observations is shown by the black filled oval in the 
lower right. 

by mass hydrogen, 27% helium, and 2% metals (Kauff- 
mann et al. 2008). Again, no uncertainties are listed 
since they are dominated by the uncertainties in mass. 

If we treat sources MMS4, MMS5, and MMS6 as a 
single group of sources and average their positions to 
determine the position of the group, the nearest neighbor 
distance for each of the 5 sources (where MMS4, MMS5, 
and MMS6 are now a single source) ranges from 0.03 — 
0.045 pc. This is comparable to the Jeans length of 0.02 
- 0.07 pc for 10 K gas at number densities n = 10 5 — 
10 6 cm -3 . However, any further interpretation is limited 
by the fact that we are only able to measure projected 
separations in the plane of the sky. 

4.2. CO 

Figure [2] displays channel maps of the SMA 12 CO J 
= 2 — 1 data at a velocity resolution of 1 km s -1 . Emis- 
sion is detected between approximately —5 to 15 km s _1 . 
The lack of emission in the channel centered around 3.25 
km s _1 suggests a cloud systemic velocity around this 
value since the cloud emission is typically fully resolved 
out in interferometer observations with similar uv cover- 
age to the data presented here. However, with no detec- 
tions of dense gas tracers, including N2D + J = 3 — 2 in 
these observations, the true systemic velocity of HBC722 
and other, nearby sources in the region is not well charac- 
terized. We note that Green et al. (2011) detected 12 CO 
J = 14 — 13 at a velocity of 6.6 km s _1 in a i7er.se/ie/-HIFI 
(Heterodyne Instrument for the Far-Infrared; de Graauw 



et al. 2010) observation centered on HBC722, but such 
a high-J transition may be dominated by warm gas in 
outflows in the region and not a good tracer of the true 
cloud systemic velocity. 

In single-dish 12 CO J = 2 — 1 data obtained at the CSO, 
Green et al. (2011) detected a bipolar molecular outflow 
centered approximately at the position of HBC722, with 
blueshifted emission extending to the southeast and red- 
shifted emission extending to the northwest. However, 
with a 32.5" beam, they lacked the spatial resolution 
to definitively identify the driving source of this outflow 
and noted the possibility that none of the detected CO 
emission was related to HBC722 itself. Figure [3] shows in- 
tegrated blueshifted and redshifted SMA 12 CO J = 2 - 1 
emission contours overlaid on the 1.3 mm continuum im- 
age, integrated over the velocity intervals of —5 to 1 
km s _1 for the blueshifted emission and 8 to 14 km s _1 
for the redshifted emission. The presence of blueshifted 
emission extending to the southeast and redshifted emis- 
sion extending to the northwest is confirmed in the SMA 
data. The emission is centered on the continuum source 
MMS3, arguing that it is this source rather than HBC722 
that is driving the outflow detected by Green et al. and 
confirmed by these SMA observations. We also note the 
presence of weak redshifted emission extending to the 
northwest of MMS5 in an L-shaped morphology (see, in 
particular, the 8.25 km s _1 panel of Figure [2]) and weak 
blueshifted emission extended to the southwest of MMS2 
(most noticeable in the 0.25 km s _1 panel of Figure [2]), 
possibly indicating the presence of an outflow driven by 
MMS5. 

The NW-SE outflow driven by MMS3 shows an S- 
shaped morphology, similar to that observed in several 
other sources over a range of masses, including IRAS 
20126+4104 (Shepherd et al. 2000), L1157 (Zhang et 
al. 2000), and RN043 (Arce & Sargent 2005). This type 
of morphology is typically explained by an outflow axis 
that precesses with time. Such precession is generally 
thought to arise from either tidal interactions between 
the disk of the source driving the outflow and compan- 
ion sources or from anisotropic accretion events (e.g., 
Shepherd et al. 2000). A similar morphology, in partic- 
ular an apparent change in the axis of the red lobe from 
one oriented more north-south to one oriented more east- 
west as distance from the driving source increases, is also 
hinted at in the single-dish data presented by Green et 
al. (2011). However, as discussed in more detail below, 
these SMA observations are only recovering the dens- 
est components of a very extended emission morphology, 
leading to significant uncertainty in the true underlying 
morphology of this outflow. We conclude that there is 
tentative but unconfirmed evidence that the axis of the 
outflow driven by MMS3 is precessing. 

High-velocity outflowing gas from MMS3 is not the 
only source of CO emission in the SMA map. Inspection 
of Figure [2] clearly shows substantial 12 CO J = 2 — 1 emis- 
sion between 1 and 8 km s _1 , the lower bounds for the 
blueshifted and redshifted emission, respectively. Figure 
Uplots both 12 CO J = 2 - 1 and 13 CO J = 2 - 1 emission 
contours integrated between 1 — 8 km s _1 . Most of the 
emission is located along the NW-SE axis of the outflow 
driven by MM S3 and likely arises from a lower- velocity 
component of this outflow. However, MMS1, MMS2, and 
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Fig. 4. — Low-velocity 12 CO (left) and 13 CO (right) integrated emission overlaid on the SMA 1.3 mm continuum image of the region 
surrounding HBC722. The background images and green contours are the same as described in Figure [3] The black contours show the 
12 CO (left) and 13 CO (right) emission integrated from 1 to 8 km s — 1 . For 12 CO they start at 3 Jy beam -1 km s — 1 and increase by 3 Jy 
beam -1 km s — 1 , while for 13 CO they start at 1.2 Jy beam -1 km s — 1 and increase by 1.2 Jy beam -1 km s — 1 . The synthesized beam sizes 
and shapes of the 12 CO (left) and 13 CO (right) observations are shown by the black filled oval in the lower right. 

-[H likely only represents the densest, most compact com- 
ponents of a more extended emission morphology. This 
statement is confirmed by Figure 2 of Green et al. (2011), 
which clearly shows large-scale 12 CO J = 2 — 1 emission 
extending over MOO" (52000 AU at the distance of 520 
pc). It is further confirmed by Figure [5j which compares 
the CSO 12 CO J = 2 — 1 spectrum at the position of 
HBC722 from Green et al. (2011) to the average SMA 
12 CO J = 2 — 1 spectrum centered on HBC722, where 
the average is weighted by the CSO beam (assumed to 
be a Gaussian with a 32.5" FWHM) and is calculated in 
practice in the image plane by adding together all emis- 
sion within 32.5" of HBC722 and down- weighting based 
on the distance from HBC722 with the assumed CSO 
beam. Comparing the CSO and SMA spectra show that 
at least 50% of the emission is resolved out over essen- 
tially all velocities for which CO emission is detected, and 
more than 90% is resolved out between 1 — 4 km s -1 . 

Ultimately, we confirm the presence of at least one 
NW-SE outflow in the vicinity of HBC722, as suggested 
by Green et al. (2011). The morphology of this outflow 
suggests MMS3 rather than HBC722 as the most likely 
driving source, and suggests that this outflow may be 
precessing. There is also an extremely tentative detec- 
tion of a weak outflow driven by MMS5. We do not find 
any conclusive evidence that HBC722 itself is driving an 
outflow. Given the uncertainties in the systemic veloc- 
ity of the region and in the true morphology of the CO 
emission due to resolving out extended emission, we are 
unable to completely rule out the possibility of additional 
outflows in the region, particularly at low velocities rel- 
ative to the uncertain core systemic velocity. Indeed, 
we note that there is 12 CO and 13 CO emission spatially 
coincident with HBC722, but lack sufficient information 
to determine whether its origin lies in one or more out- 
flows or simply in ambient cloud emission. Future obser- 
vations that provide both higher spatial resolution 




Velocity (km s~ 1 ) 



Fig. 5. CSO 12 CO J = 2 - 1 spectrum centered on HBC722 
from Green et al. (2011) (gray line) and average SMA 12 CO J 
= 2—1 spectrum centered on HBC722 (black line), with the average 
weighted by the CSO beam (assumed to be a Guassian with a 32.5 /x 
FWHM). 

HBC722 are all located near the strongest emission, and 
neither the spatial resolution nor the emission morphol- 
ogy clearly indicate whether all of this emission is due 
to the outflow driven by MMS3 or if additional weak, 
low- velocity outflows are present in this region. 

Analysis of the outflow (s) in this region is further com- 
plicated by the fact that, with minimum baselines of ~ 5 
kA (corresponding to spatial scales of 25700 AU at a dis- 
tance of 520 pc) and less than 2% of the total uv point- 
ings located at projected baselines < 10 kA (correspond- 
ing to spatial scales of 12900 AU at a distance of 520 pc), 
these observations are not sensitive to extended outflow 
emission. The 12 CO J = 2 — 1 emission shown in Figures[2] 
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Fig. 6. — SMA 1.3mm continuum intensity contours overlaid on infrared and submillimeter images of the region surrounding HBC722. 
The panels show, from left-to-right and top-to-bottom, UKIDSS K s band, Spitzer 8 ^m, Spitzer 24 /im, Herschel 70 ^m, SHARC-II 350 
/im, and Herschel 500 jum images. The cross in each panel marks the position of HBC722 itself, and the thick black lines show the primary 
beam of the SMA centered on the two pointings in the mosaic. The thin black lines show the SMA 1.3 mm continuum intensity contours, 
with the solid lines starting at 3a and increasing by 2a and the dashed lines starting at —3a and decreasing by 2cr, where the la rms in the 
continuum image is 1.65 mjy beam -1 . Pointing offsets on the order of a few arcseconds are visible between the SMA continuum emission 
and Herschel and SHARC-II continuum emission; these offsets are within the pointing uncertainties of Herschel and SHARC-II and are 
much smaller than the single-dish beams, thus we consider them to be negligible. 



better sensitivity to extended emission are required to 
fully analyze the outflows in the vicinity of HBC722. 

5. DISCUSSION 
5.1. Continuum Sources 

Figure [6] displays infrared and submillimeter images 
of the HBC722 environment, including an UKIDSS 
(UKIRT [United Kingdom Infrared Telescope] Infrared 
Deep Sky Survey) Ks image, a Spitzer 8 jam image from 
Guieu et al. (2009), a Spitzer 24 /im image from Rebull 
et al. (2011), and a Herschel 70 jam image, SHARC-II 
350 jam image, and Herschel 500 jam image from Green 
et al. (2011). Overplotted in black are the SMA 1.3 mm 
continuum intensity contours and the primary beams of 
the SMA pointings. A version of this Figure without the 
SMA 1.3 mm continuum intensity contours was previ- 
ously presented by Green et al. (2011). 

Inspection of Figure [6] shows that most of the seven 
detected SMA continuum sources are associated with 
sources at other wavelengths. A detailed discussion of 
each individual continuum source and its associations 
with sources at other wavelengths is given below. Table 
[3] lists, for each source, Spitzer photometry at 3.6 — 24 
jam from Guieu et al. (2009) and Rebull et al. (2011), 
Herschel 70 and 100 jam photometry from Green et 
al. (2011), SHARC-II 350 jam photometry from Green 
et al. (2011), and SMA 1.3 mm continuum photometry 
from this work. Figure [71 plots SEDs of all seven contin- 
uum sources, including both detections and upper limits 



(see below for more details). 

Table S] presents three evolutionary indicators for each 
source calculated from the SEDs tabulated in Tableland 
plotted in Figure [T] the infrared spectral index (a), the 
bolometric temperature (T^), and the bolometric lumi- 
nosity (Lb i). As first denned by Lada & Wilking (1984) 
and Lada (1987), a is the infrared slope in log space of 
vS v vs. v and is used to classify sources into different 
evolutionary stages (see Evans et al. [2009] for a recent 
review of classification via a). In this study we calculate 
a with a linear least-squares fit to all available Spitzer 
photometry between 3.6 — 24 jam. T^i is defined as the 
temperature of a blackbody with the same flux- weighted 
mean frequency as the source SED (Myers & Ladd 1993) 
and provides an alternative classification method to a 
(Chen et al. 1995; Evans et al. 2009). We calculate both 
Tboi and L^i by using the trapezoid rule to integrate over 
the finely sampled SEDs; a detailed description of the im- 
plementation of this method and resulting errors due to 
the finite sampling of the observed SEDs is given in Ap- 
pendix B of Dunham et al. (2008). An additional error is 
introduced by the fact that the SMA 1.3 mm continuum 
data resolve out some of the emission from the extended 
core and are thus lower limits to the true flux densities 
at this wavelength, artifically steepening the far-infrared 
and submillimeter slope of the SEDs. The magnitude 
of the error introduced depends on both the amount of 
emission recovered by the SMA and the spectral shape of 
each source; a very conservative estimate that only 1% of 
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TABLE 3 

Infrared and (Sub)millimeter Photometry of Continuum Sources 



Wavelength 


A/TA/TQ1 
IVllVlol 


A/TA/TQ9 


A/TA/TQQ 
1V11V1 oo 


A/TA/TQ/1 


i\/ri\/rcc 
IVllvloO 


A/TA/TQfi 
IVllVloO 


A/TA/TQ7 
IVllVlo I 


w 


(mJy) 


(mJy) 


(mJy) 


S u (mJy) 


S„ (mJy) 


S u (mJy) 


Su (mJy) 


3.6 


<0.6 


17.4 ± 0.9 


5.61 ± 0.28 


<0.6 


<0.6 


<0.6 


13.6 ± 0.68 


4.5 


<0.6 


25.5 ± 1.3 


10.9 ± 0.55 


<0.6 


<0.6 


<0.6 


25.7 ± 1.3 


5.8 


<0.6 


35.4 ± 1.8 


14.5 ± 0.72 


<0.6 


<0.6 


<0.6 


30.2 ± 1.5 


8.0 


<0.6 


58.0 ± 2.9 


20.1 ± 1.0 


<0.6 


<0.6 


<0.6 


38.6 ± 1.9 


24 


41.1 ± 4.1 


303 ± 30 


219 ± 22 


<2.8 


<2.8 


<2.8 


76.2 ± 3.8 


70 a 


1250 ± 130 


1260 ± 130 


6730 ± 680 


<90 


369 ± 42 


<151 


179 ± 22 


100 a 


2230 ± 230 


2250 ± 230 


9530 ± 960 


<712 


1290 ± 130 


<1020 


<458 


350 b 


2700 ± 500 


<2800 


12000 ± 2400 


<4000 


10000 ± 2000 


<3000 


<720 


1300 


38.3 ± 7.7 


55.1 ± 11 


66.1 ± 13 


60.5 ± 12 


62.7 ± 13 


44.3 ± 8.9 


20.3 ± 4.1 



a All flux densities and upper limits are derived from the Herschel 70 and 100 /im maps from Green et al. (2011). The 
flux densities are calculated in 10" diameter apertures, and the upper limits are described in the text. 

b All SHARC-II 350 /im flux densities are measured in 20" diameter apertures (see Wu et al. 2007 for more details on 
SHARC-II aperture photometry and calibration). 
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Fig. 7. — Spectral Energy Distributions (SEDs) for each of the seven detected continuum sources, consisting of Spitzer 3.6 — 24 data, 
Herschel 70 and 100 /im data, SHARC-II 350 /im data, and the SMA 1.3 mm continuum data (see text for details). Detections are plotted 
as filled circles with error bars and upper limits are plotted as open triangles. Each panel is labeled with the corresponding source. 



TABLE 4 
Evolutionary Indicators 



Source 


a 


Tbol 
(K) 


Lbol 
(L©) 


MMS1 




39 


1.1 


MMS2 


0.49 


147 


1.5 


MMS3 


0.76 


52 


5.9 


MMS4 








MMS5 




15 


1.4 


MMS6 








MMS7 


-0.20 


328 


0.3 



the emission is recovered leads to underestimates in L^i 
by less than a factor of 2 for all sources, and less than 
50% for all but MMS7. Thus, while we caution that the 
calculated L^i may underestimate the true values, the 
magnitude of these underestimates are likely comparable 
to or less than the other errors discussed by Dunham et 
al. (2008). 

Finally, we note that there are three other infrared 
sources identified as Spitzer-detected Young Stellar Ob- 
jects by Guieu et al. (2009) and Rebull et al. (2011) cov- 



ered by our SMA observations, and several other infrared 
sources covered that are not identified as Young Stellar 
Objects and are thus likely background or foreground 
stars. As our focus is on providing a high spatial reso- 
lution view of the millimeter emission in the vicinity of 
HBC722 we do not discuss these sources further except 
to note that their non-detections in our SMA 1.3 mm 
continuum observations suggest similar upper limits to 
the masses of any circumstellar disks surrounding these 
objects as for HBC722 itself, which we discuss below in 

5.1.1. MMS1 

MMS1 is associated with the Spitzer infrared source 
SST J205816.56+435352.9 from Guieu et al. (2009) and 
Rebull et al. (2011). It is detected by Spitzer only at 24 
/am with the flux density listed in Table [3j The Spitzer 
3.6 — 8 jam upper limits listed are taken from Guieu et 
al. (2009); they note that their 90% completeness limits 
increase from about 0.3 mJy at 3.6 jim to 0.6 mJy at 8 
jim but do not specifically list these limits for 4.5 and 
5.8 /im, thus we conservatively take the upper limits to 
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be 0.6 mJy in all four bands. MMS1 is also associated 
with a Herschel source detected at 70 and 100 /am from 
Green et al. (2009); all other Herschel wavelengths are of 
too low spatial resolution (18 — 35" at 160 — 500 /am) to 
accurately separate the 7 continuum sources and are not 
considered in this study. The 70 and 100 /am photome- 
try presented in Table [3] is calculated with 10" diameter 
apertures chosen as the best compromise between includ- 
ing as much source flux as possible and excluding flux 
from other, nearby sources. No aperture or color cor- 
rections are applied since neither the true spatial profile 
of the emission nor the underlying spectral shape of any 
source is known. Finally, MMS1 is also associated with a 
SHARC-II 350 /am source. The photometry presented in 
Table [3] is calculated in a 20" diameter aperture following 
the method described by Wu et al. (2007). 

The SED of MMS1 resembles that of a deeply embed- 
ded protostar. We are unable to calculate a with only 
one Spitzer detection, but the upper limits are consis- 
tent with a positive a indicative of a Class 0/1 object. 
We calculate T boi = 39 K, placing MMS1 in Class 0, 
consistent with the above statements. We thus conclude 
that MMS1 is a d eeply embedded Class protostar. As 
described in §4.21 above, we do not detect clear signa- 
tures of an outflow driven by MMS1 but do note that 
there may be additional outflows in this region that are 
not fully separated spatially and kinematically by these 
data. 

5.1.2. MMS2 

MMS2 is associated with the Spitzer infrared source 
SST J205816.8+435335.6 from Guieu et al. (2009) and 
Rebull et al. (2011) and is detected at all five Spitzer 
wavelengths from 3.6 — 24 /am. It is also detected at 
70 and 100 /am with Herschel. The photometry is again 
calculated in 10" diameter apertures and is uncertain 
since the apertures do not fully capture the emission 
that extends to the northwest but also partially overlap 
with the emission from the brighter MMS3. More accu- 
rate photometry at these wavelengths will require higher- 
resolution observations. There is no clear SHARC-II 
350 /am source, but the position of MMS2 overlaps with 
bright emission from MMS3 and MMS5. We calculate 
the upper limit as the flux in one beam at the position 
of MMS2 from the other, nearby sources. 

The SED of MMS2 resembles that of a Class I proto- 
star more evolved and less deeply embedded than MMS1. 
The calculated values of a and T^i (0.49 and 147 K, re- 
spectively) both classify MMS2 as Class I, confirming 
this statement. 

5.1.3. MMS3 

MMS3 is associated with the Spitzer infrared source 
SST J205817.7+435331.1 from Guieu et al. (2009) and 
Rebull et al. (2011) and is detected at 3.6 — 24 /am with 
Spitzer. It is also detected at 70 and 100 /am with Her- 
schel and 350 /am with SHARC-II and is the brightest 
source in the region at these wavelengths. The photom- 
etry at 70, 100, and 350 /am is calculated with apertures 
and methods identical to MMS1 above. The true flux 
densities at these wavelengths are likely higher than the 
values listed in Table [3] since the apertures do not include 
all of the extended emission, but larger apertures are 



not feasible since they would overlap with other, nearby 
sources. 

As with MMS1, the SED of MMS3 resembles that of 
a deeply embedded protostar. The calculated values of 
both a (0.76; in the Class 0/1 category) and T^i (52 K; in 
the Class categ ory) are consistent with this observation. 
As discussed in 34~2l above, MMS3 is driving a NW-SE 
outflow detected both by the SMA 12 CO J = 2 - 1 ob- 
servations presented here and single-dish 12 CO J = 2 — 1 
observations presented by Green et al. (2011). The axis 
of this outflow may be precessing over time. 

5.1.4. MMS4 

MMS4 is not associated with a Spitzer infrared source 
at 3.6 — 24 /im, it is not associated with a Herschel in- 
frared source at 70 and 100 /im, and it is not associated 
with a SHARC-II 350 /am submillimeter source. The up- 
per limits for 3.6 — 8 /am are again taken from Guieu 
et al. (2009) as described above for MMS1. For 24 /am, 
we take the upper limit to be the point at which the 
source count histogram presented in Figure 9 of Rebull 
et al. (2011) turns over, which we estimate to be at a 
magnitude of 8.5 (corresponding to a flux density of 2.8 
mJy). While there is no Herschel 70 or 100 /am source or 
SHARC-II 350 /am source, there is emission at the posi- 
tion of MMS4 from the brighter nearby sources (MMS3 
and MMS5). Thus, similar to MMS2 above, we calculate 
the upper limits as the flux in one beam at the position 
of MMS4 from other, nearby sources. 

With no detections of a compact, infrared source in any 
of the Spitzer bands, and also no detections at 70 — 350 
/im, MMS4 is possibly a starless core heated only exter- 
nally and thus too faint to detect at 350 /am above the 
emission from nearby, brighter sources. The 350 /am up- 
per limit of 4 Jy is consistent with this statement since 
starless cores and cores containing very low luminosity 
protostars are typically less than 1 - 2 Jy at this wave- 
length (Wu et al. 2007). However, recent work suggests 
that detections of starless cores with current interferom- 
eters are extremely rare since starless cores are not yet 
very centrally condensed and are thus fully resolved out 
(Schnee et al. 2010; Offner et al. 2012). If MMS4 is in- 
deed a starless core, the SMA detection indicates it may 
be very evolved and close to the onset of star formation. 
Alternatively, star formation may have already begun, 
with the core harboring a very young, very low lumi- 
nosity protostar or first hydrostatic core. Confirmation 
would require either the detection of a very faint infrared 
source below the detection limits of the data considered 
here, such as the 70 /am detection of the source Per-Bolo 
58 presented by Enoch et al. (2010), or the detection 
of a molecular outflow driven by this core, such as the 
detections of outflows from cores previously believed to 
be starless presented by Chen et al. (2010), Dunham et 
al. (2011), Pineda et al. (2011), and Schnee et al. (2012). 
No such outflow is detected in our 12 CO J = 2 - 1 ob- 
servations, but this topic should be revisited with future 
observations providing higher sensitivity and higher spa- 
tial resolution. 

5.1.5. MMS5 

MMS5 is not associated with a Spitzer infrared source 
at 3.6 — 24 /im, and the upper limits listed in Table [3] are 
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determined as described above. MMS5 is associated with 
a source detected at 70 and 100 /im with Herschel and 
350 fim with SHARC-II, and the photometry at these 
wavelengths is calculated with apertures and methods 
identical to MMS1 above. The SED of MMS5 resem- 
bles that of a Class protostar too deeply embedded to 
be detected in the mid-infrared with Spitzer. With no 
such detections we are unable to calculate a, but cal- 
culate a value for T^oi (15 K) consistent with that of a 
Class source. As noted in §4.2| there is weak redshifted 
emission extending to the northwest of MMS5 and weak 
blueshifted emission to the southeast. These weak fea- 
tures may be due to an outflow driven by MMS5, but 
higher sensitivity 12 CO observations are required to con- 
firm this tentative outflow. 

5.1.6. MMS6 

Similar to MMS4, MMS6 is not associated with a 
Spitzer infrared source at 3.6 — 24 /im, it is not asso- 
ciated with a Herschel infrared source at 70 and 100 /im, 
and it is not associated with a SHARC-II 350 /im sub- 
millimeter source. The 3.6 — 24 /im upper limits listed 
in Table [3] are determined as described above. The 70, 
100, and 350 /im upper limits are calculated as the flux 
in one beam at the position of MMS6 from other, nearby 
sources. The same discussion presented above for the 
evolutionary status of MMS4 also applies for MMS6; it 
is likely either an evolved starless core close to the onset 
of star formation or a very young, very low luminosity 
protostar or first hydrostatic core. As with MMS4, we do 
not detect any evidence for an outflow driven by MMS6. 

5.1.7. MMS7 

MMS7 is associated with the Spitzer infrared source 
SST J205817.06+435316.1 from Guieu et al. (2009) and 
Rebull et al. (2011) and is detected at all five Spitzer 
wavelengths. It is also detected at 70 /im with Herschel, 
and the photometry presented in Table [3] is calculated 
with an aperture and methods identical to MMS1 above. 
It is not detected at either 100 /im with Herschel or 350 
/im with SHARC-II; upper limits are again calculated as 
the flux in one beam at the position of MMS7 from other, 
nearby sources. The SED of MMS7 resembles that of a 
young stellar object (YSO) surrounded by a circumstel- 
lar disk but no longer embedded within a dense core; the 
calculated values of a (—0.20; in the Flat Spectrum cat- 
egory) and Tboi (328 K; near the Class I/II boundary) 
confirm this assessment. 

5.2. Evolutionary Status of HBC722 

Prior to outburst HBC722 was regarded as a Class II 
T Tauri star with a spectral type of K7— M0, a mass of 
~ 0.5 — 0.6 M , a visual extinction of 3.4 magnitudes, 
an infrared spectral index of —0.77, and a bolometric 
luminosity of 0.85 L (Cohen & Kuhi 1979; Kospal et 
al. 2011; Miller et al. 2011). As noted in g3J HBC722 
itself is not detected in our SMA 1.3 mm continuum ob- 
servations down to a 3a upper limit of 5 mJy bea m -1 . 
Under the same assumptions as discussed above in §4.1| 
this corresponds to an upper limit of 0.02 M for the 
disk mass (lower if the disk is warmer than 30 K) and 
an upper limit of 3% - 4% for the ratio of disk to stellar 
mass (again, lower if the disk is actually warmer than 



30 K). In a large submillimeter survey of circumstellar 
disks around young stars, Andrews & Williams (2005) 
report a mean disk mass of 0.005 M with a large dis- 
persion of ~0.5 dex, and a median ratio of disk to stellar 
mass of 0.5%. Thus, our observations do not rule out the 
presence of a typical mass disk. 

In a typical FUor flare with an accretion rate of 10 -4 
M yr -1 and a duration of 100 yr, up to 0.01 M of 
mass can accrete from the disk onto the protostar. For 
HBC722, however, Kospal et al. (2011) calculated a burst 
accretion rate of only 10 -6 based on their measured L^i 
during the burst 12 . Our derived upper limit for the disk 
mass of 0.02 M is thus consistent with providing a suf- 
ficient mass reservoir to support the observed outburst 
in HBC722 unless the accretion rate is several orders of 
magnitude higher than estimated by Kospal et al. and/or 
the burst duration is much longer than the typical 100 
yr. Even if one of these cases were true, there is possi- 
bly gas remaining in the vicinity of HBC722 that could 
still accrete onto the star+disk system and power the 
burst, as suggested by 13 CO J = 2 — 1 and Herschel 
far-infrared continuum emission spatially coincident with 
HBC722. Further constraints on the amount of circum- 
stellar mass available to accrete onto HBC722 and the 
likelihood of sufficient mass remaining to power further 
bursts beyond the current one require deeper millime- 
ter continuum data probing to lower disk masses and 
higher-resolution far-infrared and submillimeter contin- 
uum data observed during the burst to better determine 
the burst luminosity and implied accretion rate. For the 
former, the high sensitivity of full-science ALMA opera- 
tions will be the ideal facility despite the high declination 
of HBC722 (+44°) since, according to the ALMA sensi- 
tivity calculator 13 , even a short, 30 minute track with 
50 antennas will improve the mass sensitivity by a factor 
of 100. For the latter, the upcoming 25-m submillimeter 
telescope CCAT will be of particular value assuming the 
burst is still in progress when CCAT begins science op- 
erations (currently expected in 2015 - 2017; Radford et 
al. 2009; Sebring 2010). 

6. SUMMARY 

In this paper we have presented 230 GHz Submillime- 
ter Array continuum and molecular line observations of 
the newly discovered FUor candidate HBC722. We sum- 
marize our main results as follows. 

1. Seven 1.3 mm continuum sources are detected in 
the vicinity of HBC722; none are HBC722 itself. 
We compile infrared and submillimeter continuum 
photometry of each source from previous studies 
and conclude that three are Class embedded pro- 
tostars, one is a Class I embedded protostar, one 
is a Class I/II transition object, and two are either 

12 This measurement is rather uncertain due to the variability of 
the outburst brightness since the initial flare in 2010 (see 32} an d 
the fact that Kospal et al. lacked photometry during the burst at 
A > 10 /im, although the latter point is mitigated by the fact that 
the Herschel 70 fim flux density of 412 mJy measured during the 
burst and reported by Green et al. (2011) is generally consistent 
with the assumptions made by Kospal et al. to extrapolate beyond 
10 /im. 

13 Available at https://almascience.nrao.edu/call-for- 
proposals / sensitivity-calculator 
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starless cores or very young, very low luminosity 
protostars or first hydrostatic cores. 

2. A northwest-southeast outflow is detected in the 
12 CO J = 2 — 1 observations. This outflow is cen- 
tered on and thus likely driven by MMS3, one of the 
Class sources detected in the 1.3 mm continuum 
data, and its axis may be precessing. This outflow 
detection is consistent with a similar outflow de- 
tected in low-resolution, single-dish 12 CO J = 2 — 1 
observations presented by Green et al. (2011). Our 
higher spatial resolution confirms that HBC722 is 
not the driving source. 

3. There is no conclusive evidence that HBC722 it- 
self is driving an outflow, although we caution that 
higher spatial resolution, better sensitivity to ex- 
tended emission, and better determinations of the 
systemic velocities of the sources in the vicinity of 
HBC722 are needed to fully evaluate the kinemat- 
ics of the 12 CO J = 2 — 1 gas in this region. 

4. The non-detection of HBC722 in the 1.3 mm con- 
tinuum observations suggests an upper limit of 0.02 
M for the mass of the circumstellar disk, consis- 
tent with typical T Tauri disks. This upper limit is 
consistent with a disk that provides sufficient mass 
to power the burst. Future observations are needed 
to further study the actual amount of circumstel- 
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lar mass available to accrete onto HBC722 and the 
likelihood of sufficient mass remaining to power ad- 
ditional bursts beyond the current one. 

We have noted in the text several future observations 
that are needed in order to better disentangle the mil- 
limeter emission in this complicated environment and 
better determine the properties and evolutionary status 
of HBC722. 
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